A formally exact field theory for classical systems at equilibrium. 



D. di Caprio * and J. P. Badiali 

Laboratoire d'Electrochimie et de Chimie Analytique, UMR 7575, 
UPMC Univ Pans 06, CNRS, ENSCP, 
Universite Paris 6, Bat. F74, P.P. 39 
4, Place Jussieu, 75230 Paris Cedex 05, France 

We propose a formally exact statistical field theory for describing classical fluids with ingredients 
similar to those introduced in quantum field theory. We consider the following essential and related 
problems : i) how to find the correct field functional (Hamiltonian) which determines the partition 
function, it) how to introduce in a field theory the equivalent of the indiscernibility of particles, 
Hi) how to test the validity of this approach. We can use a simple Hamiltonian in which a local 
functional transposes, in terms of fields, the equivalent of the indiscernibility of particles. The 
diagrammatic expansion and the renormalization of this term is presented. This corresponds to a 
non standard problem in Feynman expansion and requires a careful investigation. Then a non-local 
term associated with an interaction pair potential is introduced in the Hamiltonian. It has been 
shown that there exists a mapping between this approach and the standard statistical mechanics 
given in terms of Mayer function expansion. We show on three properties (the chemical potential, 
the so-called contact theorem and the interfacial properties) that in the field theory the correlations 
are shifted on non usual quantities. Some perspectives of the theory are given. 



I. INTRODUCTION 



In various domains of physics a description in terms of fields is frequently used. Hydrodynamics represents a 
first example in which some fields (densities, velocities, ...) are introduced for describing the properties of a coarse 
grained entity - the so-called fluid particle. Later, field theory (FT) has been used as a simple and intuitive tool 
to predict the behaviour of complex systems in the domain of soft matter physics [H, 0, HI, El HI- These FT are 
essentially phenomenological and rely on Hamiltonian functionals introduced in an ad hoc manner. They focus on a 
mesoscopic scale description and they are based on a more or less explicit coarse graining procedure. In this context, 
Hamiltonians are introduced to describe large classes of phenomena having similar properties though different in their 
microscopic details. This further suggests another type of problems where a FT is extensively used i.e. the description 
of critical phenomena. Here also the FT is based on the assumption that a detailed microscopic knowledge of the 
system is not relevant to describe its universality class [H, 0, S, Q- And well suited approximations to describe 
systems with long range correlations or interactions are introduced. In relation, field theory is also used to describe 
systems with the long ranged Coulomb interactions. In this case, FT is constructed using the Hubbard-Stratonovich 
transform of the standard partition function [l(| EH, EH . Better known as the sine-Gordon transform in the case of 
the Coulomb potential, it has given rise to considerable literature [H, E3, EH EH E3, El, Ell US H3 • These approaches 
give an exact description of the systems properties on a microscopic level. In this respect, they are distinct from 
the soft matter like descriptions based on a coarse graining procedure. The sine-Gordon approaches introduce an 
auxiliary field and intricate couplings between fluctuating fields. In our opinion, this auxiliary field is essentially a 
mathematical tool, difficult to associate with any physical quantity. As a result finding meaningful approximations is 
rather counterintuitive and the application of such approaches requires that one be rather cautious [22j |. 
In contrast to these approaches, our main goal is to show that it is possible to write a FT directly in terms of fields 
using methods similar to those used in quantum field theory (QFT). Namely, we show that it is possible to build the 
theory around a field, which is a real quantity having a simple physical meaning. Moreover we will show that our 
FT construction is not only simple and intuitive but also leads to a complete description at microscopic level. In this 
paper, we consider systems at equilibrium. 

The paper is organized as follows. In the next Section we present the main requirements which a FT must verify. In 
Section Hill we give the Hamiltonian on which the FT is based : it contains two terms of different nature. This leads 
to investigate the Feynman expansion of a purely local Hamiltonian with an infinite number of coupling constants. 
This is developed in Section llVI where some important specific aspects of the expansion are shown. In Section fVl we 
calculate the partition function in the presence of an interaction pair potential. In the next Section we establish an 
exact mapping between our FT and the standard statistical mechanics given in terms of Mayer expansion [23| . In 
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Section lVHI we illustrate on several examples how our approach may lead to new aspects in statistical physics. Finally 
in Section [Villi we give some conclusions and perspectives. 



Our main assumption is that the partition function 3[0] of a classical system can be described exactly via a 
functional integral defined according to 



in which is a field, H[cj>] a functional of this field which we call Hamiltonian and /3 — I/(/cbT) is the inverse of the 
temperature. 

To use ([l} we must solve several problems. We have to define <fi and also find an explicit form for H [</)]. It is intuitive to 
choose for <p a real quantity as the density of matter p for instance. This choice represents a fundamental difference from 
the Hubbard-Stratonovitch type approaches in which two fields are used, one being a complex quantity. In comparison 
with the standard description of the liquid state, where the configuration space spans all possible distributions of the 
particles, here p is a function defined everywhere in space. As a consequence the number of degrees of freedom in ([T]) is 
related to the space discretization required to calculate the functional integral as opposed to the number of particles. 
Then we have to solve a new problem, of how to transpose in FT a property mimicking the indisccrnibility of particles. 
To relate the FT and the usual physics we assume that the average of the field, < p >, corresponds to the actual 
density of particles noted p 1 . 

In addition to the indiscernibility of particles, the so-called classical statistical mechanics contains the volume of 
the elementary cell ApAx = K or, at least, after integration over momenta in the case of systems of at equilibrium, 
the thermal de Broglie wavelength A. Thus we also have to find how such quantities associated with particles appear 



In contrast to the Hubbard-Stratonovitch type descriptions which are finite as they represent a rigorous mathe- 
matical transformation of a finite quantity, we expect that a field theory, as is often the case, will contain infinities 
associated with the short scale spatial discretization of the functional integral. This will require the introduction of 
a renormalization procedure. 

Finally, to be able to assert that the FT is also an exact representation, we have to show that there exists a rigorous 
mapping between the FT and the standard statistical mechanics of dense systems. Hereafter we turn our attention 
to all these questions. 



1 Quantities associated with the thermodynamics as opposed to those calculated from the FT will generally be indicated with a tilde. 



II. REQUIREMENTS FOR A FIELD THEORY 




(1) 



in a FT. 
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III. DEFINING THE HAMILTONIAN 



To build H[p] we follow an approach inspired by the methods developped in QFT where instead of H[p] the 
lagrangian L[<p] is considered. To find the latter, we select a functional and check that the mean field approximation 
of the theory reproduces a well known result; for instance, the Maxwell equations in quantum electrodynamics. In 
this case, L[</>] can be considered a good choice for elaborating the complete theory including the fluctuations via the 
functional integral. 

In statistical thermodynamics, it is only for systems without interactions, ideal systems, that we know an exact 
and general result and have the explicit expression of the partition function, So- We then consider such systems and 
require that the Hamiltonian Ho[p] reproduce the thermodynamic partition function So in a mean field estimation 
of ([1]). However, since So is a cornerstone in classical statistical mechanics and contains fundamental physics related 
to A and the indiscernibility of particles, we further require that Hq [p] reproduce So exactly, i. e. also beyond the 
mean field approximation. These fundamental aspects will then be correctly accounted for in the FT for all systems 
including those with interactions. We shall now discuss the Hamiltonian. 

Having chosen the field p so that its average corresponds to the density of particles p, it is evident to fix the chemical 
potential, p, and choose for So the grand canonical partition function. In this case we have the exact thermodynamic 
results 

lnE = [3pV = pV (2) 

and 

(3p = ln(pA 3 ) (3) 

where the exact density p is uniformly distributed in space. A simple Hamiltonian Ho[p] that reproduces ([3]) in a 
mean field approximation of ([1]) is 



(3H [p] = J dr (p(r) [in (p(r)A 3 ) - l] - j3pp{v) } . (4) 

As for the lagrangian L[<j)] in QFT we cannot claim that ifo[p] 1S unique. However we see that the part 

F[p}= J drp(r)[ln(p(r)A 3 )-l] (5) 

of Ho[p] represents, in the mean field approximation, the free energy of an ideal system i.e. the kinetic energy and 
the entropy. In Section l"VII[ we compare F[p] with the DFT (density functional theory) [13, HI, HH| where a similar 
term appears. 

The requirement that Hq[p\ gives the exact result entails a more careful analysis of the functional integral So[p]. 
In order to calculate practically |T]), we have to introduce in the r.h.s. a lattice with a spacing a. The result will then 
depend on this parameter. In the following, our intention is to find the conditions to obtain the exact thermodynamic 
result whatever the value of a. The discrete form of [3H [p] is 

V/a 3 V/a 3 

f3H Q [p] = ]T p(r,)a 3 [in (p(ri)A 3 ) - l] - £ p(r,)a 3 (6) 



and the partition function becomes 



r V/a 3 

So[p] - / H d[p(^)a 3 } e-^M (7) 



where we have used in the measure the dimensionless quantity p(vi)a 3 . Due to the local character of Ho[p], the 
calculation of So is a product of usual integrals, such as 

d[p(r 2 )a 3 ] exp {-p^a 3 [in (p(r,)A 3 ) - l] + f3pp{T l )a 3 } . (8) 
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Beyond the saddle point, we have 

lnS [p] = pV 



1 

pa 



~ 3 4>[pa 3 } 



(9) 



where the function ip given in [X] represents the correction to the exact thermodynamic result pV given in As 
expected, the correcting term becomes negligible when pa 3 is large. However, the discretization has introduced a 
cumbersome dependence on the lattice spacing a, which we would like to dispose of keeping only physically meaningful 
terms. Before discussing this point, we generalize this result for an external potential. 

Having to calculate local integrals the previous result for the partition function can be easily extended by changing 
(3p into (3p — V ext (ri) where the external potential is in temperature reduced units. Instead of ([9]) we now have 

V/a 3 V/a 3 

lnSo[p, V ext ] = /5a 3 e~ yC " (r0 + £ ^{pa 3 e~ v ^^} (10) 

i i 

where the last term in the r.h.s. of ([10]) is given in[XJ If V ext (ri) 1, the corrective term is still negligeable when pa 3 
is large. Equation (fT0|) is correct as long as V ext {vi) does not vary rapidly on the distance a which is already large 
in comparison with the mean distance between particles (~ p -1 / 3 ). This condition is a restriction on the validity of 

m- 

It is possible to release such constraint and generalize eq. (|10p to any external potential by noticing that all physical 
terms for this system have a well defined dependence on the lattice spacing, distinct from the corrections associated 
with ip. We now take into account the following quantity 

V/a 3 

lnZ«{p,V ext ] = lnE [V ext }-Y, ^ e~ v ' mt ^] (11) 

i 

V/a 3 

= $> 3 e-^). (12) 

i 

The renormalized quantity is now equal to its value at the saddle point whatever the value of a. From an operational 
point of view this result must be understood as follows : exp{— f3Ho[p]} is a formal expression. It represents an expan- 
sion of the exponential around its saddle point value and in this expansion terms corresponding to tp{pa 3 e~ v ^ ri '] 
are discarded in order to obtain the physical quantities. 

In the limit a — > 0, the renormalized grand potential is now a finite quantity and its value 

ln~«{p,V ext }= f~pe- v ^dv, (13) 



corresponds to the standard statistical mechanics result valid for any external potential which is independent from a. 
Note that the change of limit due to the presence of a potential is a standard problem in statistical mechanics as 
shown in [27[ • For an ideal system the so called classical statistical mechanics is obtained in the limit h — > whatever 
the value of A. However, in the presence of an interaction potential, an extra limit A — > must be taken in order to 
keep all information about the interaction potential. 



From the results obtained in this Section we assume that the total Hamiltonian will be in the form of 

PH[p] = PH [p] + l -J dTdv' p{T)Pv{v - r')p(r') (14) 

where at this stage there is a non-local term due to the presence of the interaction pair potential v(r — r'). 

Hereafter our main goal is to give an operational meaning to (fT4")) as we have already done in the case of (3Ho[p\. 
To calculate the partition function, we need to expand exp{— (3H [p]}. In QFT the calculation of similar quantities is 
done by introducing a gaussian propagator and performing the so-called loop expansion. Here Hq[p] is purely local 
and it is not traditional to give a Feynman expansion for such a term. We have to find a formal propagator and a loop 
expansion associated with Ho[p] in order to be able to treat the local and non local part of H [p] on the same footing. 
In the next Section we shall investigate the properties of H [p] and will show that this expansion is also fundamental 
when an interaction potential is present. 
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IV. FEYNMAN EXPANSION OF H [p\. 

By using the fact that Hq[p] is local we have obtained a first expansion cqs. (fT0|) and (|A5|) of SolV 61 *]. In parallel a 
second expansion can be performed with Feynman graphs. Since both expansions are expressed in terms of pa 3 and 
_yext we ma y identify term by term the expansion coefficients and by this method, as we shall see, solve complex 
problems of combinatory. 

To have a simple Feynman expansion, we choose a constant field p z as a reference state. Our choice - obviously not 
unique - is to take p z = e^ M /A 3 corresponding to the activity [28] 2 . This choice combines two interesting points : i) 
it gives an expansion in terms of activity which will be useful when comparing in Section IVII our results to the Mayer 
expansion and ii) it leads to a very simple propagator. 



A. Gaussian propagator, perturbative expansion. 



Hereafter we write the field as p(r) — p z + Sp(r) and we have 

/3H [p,V ext ] = (3H { 0) [p z ,V ext }+(3H { 2) [Sp}+f3SH[Sp,V ext } (15) 



where the first term is constant 



V/a 3 

(3H^ [ Pz , V ext ] = ~ Pz V + E V ext (n)p z a 3 . (16) 



The second term is quadratic 

V/a 3 



= -i-j E Spiny % Sp^a 3 (17) 

and the Kronecker £y plays formally the role of an interaction. Following the terminology of the QFT we call the 
quadratic term propagator. The remaining term represents the coupling part of the Hamiltonian given by 

v/a3 v/ai 00 r-iW/-2V /i \ 

[3SH[s P ,v ext ] = E v ext (n)s P (n)a 3 + E E \ Pz J)v-i) {nVpfryi 1 ) ■ ( 18 ) 

It shows the specificity of the present FT, with an infinity of coupling terms whose coefficients depend on a numerical 
factor and the parameter p z a? . 

In order to achieve a diagrammatic expansion we rewrite the partition function according to 



Zo[V ext ,J] = J Y[d[p(n)a 3 } ex v\~ a +E J ( r ')^) fl3 



where Ho[p,V ext ] includes the external potential, J is a generating field and a is a parameter, formally equal to 

in QFT 0, §. We can perform the 
t formally as 

PH^[ Pz ,V^}\ , r— 



genere 

1 , which is useful to organize the loop expansion in QFT [3, Q . We can perform the functional integral, using the 
gaussian integrals [f| @, ,8, 9] and express the result formally as 



5 [^,J] = ex P |- ^° ^^ (V^pTa 
exp [Yr VeXt ] jexpj^E J(n)kyl{* 



(20) 



2 p and p z are identical in the present case of the system without interactions. 
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The second line introduces the operator obtained by replacing the field Sp(r)a 3 with the derivation operator 5/5 J(r) 

This operator is applied to the last term in the r.h.s. of ((20)) which is gaussian [1, 0, H, [9[ ■ The calculation is performed 
expanding the operator exp {— (f3/a)SH {5/5 J, V ext }}. Taking J = at the end of the calculation, we select terms 
with pairs of derivatives acting on the same quadratic form, which corresponds to the well known Wick theorem 
@) S B ■ Note that going from the density field representation of 2 to the generating functional representation, 
we invert the kernel of the quadratic form. In the present case, this is a simple operation which consists in taking the 
inverse coefficient. 



B. Diagrammatic representation 

The diagrammatic representation of the theory is organized [1, 0, H, Q around vertices representing couplings of 
the field obtained from the expansion of exp{— [35H/a} and lines joining the vertices representing the propagator. 
The symbols used to draw these elements are shown in figure [TJ The propagator is represented by a curly line. The 



n > 3 coupling term external weight propagator 

FIG. 1: Diagrammatic elements for the graph representation. 

coupling terms will be denoted by a black circle, where three or more propagators can be joined, the precise number 
is understood from the number n of lines joined to it. An important feature in knowing explicitly the Hamiltonian 
functional is that, as opposed to the phenomcnological FT, we can workout precisely all coefficients for the couplings. 
Hence, besides the standard l/n\, the coefficient is (— l)( n_:L )(n — 2)!/( / o z a 3 )"~ 1 . The case of the vertex, with only 
one line attached to it, is drawn by a crossed circle and is associated with the external weight — V ext (r). Depending 
on whether we use the generating functional representation or not, a coupling term may represent [5p(r)a 3 ] n or 
[5/5J(r)] n . Furthermore, otherwise specified, we shall take into account connected graphs related to the logarithm of 
the partition function. 

Let us now define some topological elements. The external branches are the one body coupling constants together 
with the only propagator which can be attached to it. Internal lines are propagators which are not in external 
branches. 

The graph in figure[2]is an example of a diagram. The points 7, 8, 9 are vertices, the points 1, 6 are the external 




FIG. 2: Example of a diagram, 
weights and there are 6 external branches and two internal lines. 



7 



C. Topology. 

1. Dimensional analysis. 

The diagrammatic representation of In So [p] leads to an infinity of graphs that we can classify, as common practice 
in FT, by using a dimensional analysis in terms of the parameter a 6, 7, 8, 9]. This corresponds to the loop expansion. 
A diagram with L loops is dimensionally associated with a L_1 . For instance, the diagram in Figure [2 which is a 
tree diagram (L = 0), is indeed proportional to 1/a. In the following, our purpose is to show that the dimensional 
analysis in terms of the formal parameter a can be associated with a physical parameter of the system. 
The standard analysis allows to relating the number of elements in a graph (lines, vertices) to the number of loops 
[1, 0, 0, 11 according to 

L-l = I + E-m (22) 

where I is the number of internal lines, E the number of external lines, and m the number of vertices. These include 
also the one point vertices. The latter associated with the external potential set the power in —V ext . The r.h.s. of 
this relation shows that the power in a of the graph corresponds, in agreement with the expression of the partition 
function eq. ([20]) . to a factor 1/a for each of the m vertices, a^ I+E ^ for the internal and external lines. It is tempting 
to consider the quantity p z a 3 instead of a which appears in the calculation in a similar way. Here, we also have to 
account for the power of this term in each coupling term, eq. (|18p . Let rm be this power for each of the m coupling 
vertices. Considering that each line is attached to two vertices, we have 

m 

rm = 21 + 2E. (23) 

i=l 

Using eq. ^ 

m 

^Tmi = 2(L + m - 1) (24) 
i=i 

which is also 

m 

-J2(rrn-1) + I + E = -(L-1) (25) 
i=l 

On the left hand side we recognize the contribution in powers of p z a 3 in the graph : each vertex contributes 
l/(p z a 3 ) mi ~ 1 , and there are I + E lines each contributing p z a 3 . Thus we obtain a relation between the overall 
power of p z a 3 of the graph and the number of loops. The role of the parameter for p z a 3 is then equivalent to 
that of 1/a. In the following, we no longer introduce the factor a, as its role is redundant. All graphs can be 
computed exactly and we shall avoid explicit indexing of the points in the expressions, as finally all points are the same 
and we shall only discuss combinatory. The value of a graph is a numerical coefficient, a power of p z a 3 and of — V ext (r). 



2. Tree graphs. 

First we take the case L = 0, which corresponds to tree graphs. From dimensional analysis, all tree graphs with n 
external branches are proportional to p z a 3 [— V ext (Vi)] n and the value of their sum can be written as 

Pz a 3 c n [-V^inY (26) 

where c„ is a combinatorial coefficient. The value can be obtained by equating this expression, linear in p z a 3 with 
the corresponding term in eq. (jlip for each point r,. Thus 

Pz a 3 (f; c n [-V ext (r t )]A = p z a 3 e- vext ^. (27) 

\n=0 / 

Order by order in powers of V ext , this equation sets c„ = 1/n! for any n > 3 and we can generalize the notion of 
trees to all n. Indeed, one can verify that the cases n = and n = 1 relate to the expression of exp {— f3H[p z a 3 \\ in 
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eq. (|15p and give cq — 1 and c\ = 1 moreover the calculation of the quadratic term in V ext gives = 1/2. Now we 
know the combinatory for the rt-tree graphs. The result is extremely simple. 

Of course c n can be also calculated directly by performing the sum of graphs, such a direct calculation shows that 
the rather simple and intuitive value of c„ results in fact from the combination of different graphs. 



3. Loop graphs. 

Let us now consider the class of connected diagrams which have at least one loop (L > 1) and n external branches 
that we refer to as n L -loop graphs. 

For L > 1, the dimensional analysis states that a given L corresponds to a power of l/(p z a 3 ). We consider graphs 
with L > 1 loops and n external branches. For given values of L and n, the dimensional analysis for the sum of all 
such graphs leads to the expression 

dh ° n ' L < -V ext {Yi)] n (28) 



where dr, are the coefficients of the expansion of ip given in[^J and c' n L is a combinatory coefficient. 
The case L = 1 is specific, for ti = 0we have 

iln(2 W 3 ) (29) 

and for n > 1 

fkl[_V^(r 4 )]". (30) 
n\ 

The contributions for any n in eq. (|28|) - (|30[) can be obtained from term by term identification with the function ip 
in eq. (|A5|) for all powers of p z a 3 and —V ext . For L = 1, the comparison gives c[ 1 = 1/2 and for n ^ 1 and c' n x = 0. 
For L > 1, we must have L = (1 — L) n . Clearly, the c' n L can be also calculated by performing the sum of the 
corresponding graphs. 



D. Ideal system vertex functions. 



In the following, we define an important object in the diagrammatic expansion. For n ^ 2, we define the n T -vertex 
functions as the functions obtained from n-tree graphs by erasing the n external branches 3 . The value of the sum of 
all graphs contributing to a n T -vertex function is 



(p z a*) 6J(r) 



(31) 



where l/(p z a 3 ) n derives from the fact that we have erased from the tree graph n external propagators and [5/6J(r)] n 
refers to the n points where this vertex function can be combined to the rest of the graph. The combinatorial coefficient 
is that of the corresponding n-tree graph. The generalisation for the case n = 2 will be given later. The general 
expression is applicable in this case also. The expression of these tree vertices constitutes an important result of this 
paper. It states that despite the variety of graphs contributing to an n T -vertex, all occurs as if we have a standard 
coupling of the field at a given point, with a coefficient which besides the standard 1/n!, is simply 1. 

Starting from graphs with any number n of external branches and loops L, we define the n L -vertex functions, 
obtained similarly to the tree vertex functions, by removing the n external branches. For L = 1 there is a single non 
zero term for n = 1 



(32) 



3 Note that these vertex functions are not the 1-particle irreducible functions of the field theory associated with a Legendre transform 
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the other terms for u / 1 are zero. And for a given L > 1 and n, the value is given by 



(1 - L) n 



1 



(PzO- 3 ) 



3YL-1 n \ 



1 



(p 2 a3) <W(r) 



(33) 



where the coefficient is that of the n^-loop graphs. The term in square brackets is again simply related to the fact 
that we have removed the n external branches and created the corresponding attaching points. 



E. Renormalization. 



In the previous analysis, we have associated topological properties of n-tree and n L -loop graphs to given powers of 
p z a 3 and of —V ext . This has been done in order to relate further this topological analysis with the analytic expression 
of the generating functional eq. (fTUl) . The sum of tree graphs corresponds to the first term in this equation, whereas 
graphs with at least one loop are part of the second term. As mentionned in Section IIII1 expression (jlOp depends 
on the lattice spacing whereas the interest of the renormalized partition function eq. (jlip is that it has a finite limit 
independent of a for vanishing lattice spacing. 

Here, in order to free ourselves from the lattice spacing and obtain the renormalized partition function, we define the 
following renormalization procedure which consists in subtracting all graphs with at least one loop. This is equivalent 
to susbstracting the term corresponding to the function ip in the analytic expression of the partition function, eq. (|10j) . 
This procedure gives a meaning to the formal expression exp{— f3Ho[p]} by giving an operational description in terms 
of diagrams. Note that after renormalization, we no longer, strictly speaking, consider Hq and thus this functional 
should not be directly compared with other formalisms where it appears. From this procedure, we now have a 
diagrammatic expansion of the renormalized partition function which corresponds to the exact result for an ideal 
system and which can be used for any value of a in particular in the limit of vanishing a, which we discuss in the next 
Section. 

In the following, we shall study the system with interactions and show that the same graphs as discussed in this 
Section appear. We will see that for the reason of locality the renormalization described here can be applied in this 
context and that we can obtain a well behaved theory also for the system with interactions. 

The present discussion may appear like a cumbersome way of treating the simple ideal system. However, the crucial 
point is to understand how the counting properties for the particles transpose to the FT. In the following, the main 
tools introduced in this Section will be used to analyse the case of a system with interactions, as we are now able to 
expand in the same way both local and non local terms in H[p] using the Feynman expansion. 
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V. FEYNMAN EXPANSION OF THE FULL HAMILTONIAN 



Hereafter we study the generating functional 

V/a 3 ( V/a 3 ~\ 

E[J] = J J[ d[p( ri )a 3 ] exp j -PH[p] + £ J(r,)p(r l )a 3 | (34) 

in which H[p] is given in (12). Expanding the field around the activity p z , we obtain 

(3H[p\ = f3H^[ Pz ,V ext }+ (3H^[Sp] + [3SH[S P ,V ext ]. (35) 

The first contribution is 

0H^[p z ,V ext ] = - Pz V+^p z Vv +p z V ext (36) 

where we define vo — (3p z Y^j^u v ( r ij) a3 an d Vo xt = V ext (ri)a 3 . In the following, we assume that we have 

substracted the self energy and that the interaction potential cannot be taken at the same point, although to simplify 
the notation we do not explicitly indicate it. The quadratic Hamiltonian is 

V/a 3 

0HW[8p] = 3 V Sp(ri)a 3 [% + 0p z a 3 v( rij )} Sp(r,)a 3 (37) 

hJ 

As noted earlier, the Kronecker Sij will be treated as an interaction. The coupling Hamiltonian is given by 

V/a 3 V/a 3 

pSH[Sp,V ext ] = £ Spavin) +J2v ext (r i )6p(r i )a 3 

i i 

where v(r) = (3p z a 3 v(v). We point out that this coupling Hamiltonian is essentially the same as for Ho[p] with 
the exception of a linear term which includes the interaction potential. Therefore the topology of the diagrammatic 
expansion will be similar to the expansion for Hq [p] . The main modifications are in the existence of a new contribution 
to the propagator and to the one body term. We then have for the generating functional 

E[J] = C xp{-/3H^[p z ; J]} ( v /2^a 3 ) V ' /a (39) 
f 3 y / a3 



exp 



-j3SH 



hi 



where, like in Section HVl we have substituted Spa 3 by the S/5J and the notation [. . .] 1 indicates the inverse. The 
latter can be expanded according to 

V/a 3 

[Sij + v{ri,j )] _1 = 5^ - v(nj) + ^ v(nk)v(r k j) + .... (40) 

fc 

In this expression, the Kronecker S^ is its own inverse and the rest represents a sum of terms of alternate signs 
constituted with chains of single potentials. The diagrammatic representation of this equation is given in figure [31 
where the full propagator appears on the l.h.s. while on the r.h.s. the curly line is the Kronecker Sij and the lines 
represent a single interaction potential — u(r,j). We can thus generalize the notion of the tree vertex function of 
Section TlVDI to the two body coupling term associated with a weight l/(2p z a 3 ). The diagrammatic expansion will 
be the same as the one given in the previous section, except that the full double line replaces the curly line and that 
the external weight has two contributions shown in figure [4j 
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+ ... 



FIG. 3: Diagrammatic representation of the decomposition eq. (1401) . 



® 



-v(r 




FIG. 4: Diagrammatic notation for the weights at the end of the external branches. 



A. Topological reduction : ideal system vertex functions. 



Hereafter we expand the propagator according to the decomposition shown in figure [3] The purpose is to apply 
a topological resummation of the theory in terms of the vertex functions introduced in Section HVDl These vertex 
functions include at least two attaching points. The case of the one body coupling term will be detailed separately. 

On the graph given in figure El we present an example of this expansion, where the diagram on the right represents 
a possible decomposition of the total propagator of the original graph on the left. We have omitted the labels and 
arbitrarily chosen one of the external weights — V ext . On the right, for simplicity, we have chosen only the contribution 



to the chain of interactions corresponding to a single interaction. These aspects are irrelevant to the present discussion. 
Given the local nature of the ideal system couplings and propagators, it is interesting to isolate in the diagrams the 
local parts which are indicated inside the dotted frame on the figure. Their contribution to the graph is a numerical 
coefficient as they are independent on the rest of the graph. 

We then consider graphs with the same backbone structure in terms of the interaction potentials but with a different 
local part like, for instance, in figure [6l The sum of all such local diagrams can be performed using the n T -vertex 
functions as defined in Section FlVDI In the present case, it requires the 5 T -vertex and 7 T -vertex functions derived from 
the 5 T -tree and 7 T -tree. The resummation into vertex functions is equivalent to a topological reduction. Noting the 
ideal system vertices introduced in IIV Dl by black squares, the graph is now represented by figure [7] Clearly, the new 
graph corresponds to various topologically different graphs of the original expansion in terms of the total propagator. 
The factor associated with each ideal system vertex function is simply according to Section HVDI : l/[n!(p z a 3 )( n ~ 1 ' ) ]. 
We hereafter detail the special case of the one body coupling which according to figure H] can have different weights : 
the interaction potential or the external potential, i) First we discuss the case when the external branch has a weight 
v. This weight can be attached on the ideal system propagator or on a chain of one or more potentials. These two 
cases are complementary to allow for any number of interactions in the chain 4 , ii) Then we consider the external 
branch associated with the external potential. The case when there is a single interaction potential on which we 
attach the external potential is specific. It is shown in figure [8] where one may verify that all cases with any number 
of — V ext are represented. For the case of a single interaction potential and all other cases, it is straightforward to 
find that any coupling term can be decorated with a factor e~ v 

The topologigal reduction, presented here on a tree graph using the ideal system tree vertex, can be generalized 
to graphs which include loops for the ideal system. The separation into non local and local parts related with the 
ideal propagator can be performed in an identical way. We then also need to introduce the ideal system loop vertices. 





FIG. 5: Example of a tree diagram with one possible decomposition. 



4 We recall that from the decomposition in figure [3] there cannot be an ideal and an interaction propagator in series. 
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FIG. 6: Diagram which has the same structure for the interactions as the diagram in figure [S] but a different topology for the 
ideal propagator. 




FIG. 7: Representation of the sum of all diagrams with the same structure for interactions as the diagram in figure [5] using the 
5 T -vertex and 7 T -vertex functions. Amongst these we have the diagrams of figure [5] and figure [6] 



Having performed the topological reduction of all local parts using either tree vertices or loop vertices, we note that 
graph with an identical structure may appear once with a tree vertex and once with a loop vertex. For a given position 
of all other vertices, we can combine the tree vertex and a loop vertex as they are taken at the same point knowing 
that the rest of the graph is identical. The sum of these two vertices corresponds to the two terms in the ideal system 
partition function eq. (| 10|) . At this stage, we can introduce the renormalization presented in Section llV El which 
corresponds to substracting the loop vertices. As a result we have vertices which are well behaved in the vanishing 
lattice space limit. 

As the potential couples distinct points, note that there cannot be a loop consisting of a single interaction potential. 
More general loops which may include chains of interaction potential are not concerned by the topological reduction 
associated with the ideal system and remain unchanged. 




two or more external potentials 



FIG. 8: Resummation of external branches with external potential on the one body coupling terms. The gray circle may 
represent any diagram and the black square a three or higher vertex function. 
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B. Diagrammatic definition of the partition function. 

The result of this Section is that the logarithm of the partition function is given by all possible connected graphs 
made of non labelled ideal system vertices i.e. n T -vertex functions [n > 2) and internal lines corresponding to a single 
potential. The coefficients of the vertices are those of the n T -vertices given in Section IIVDI The external branches 
are either the ideal system or a single potential propagator. At the end of the external branches we find the labelled 
weights shown in figure [U 
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VI. MAPPING BETWEEN THE FT AND THE MAYER EXPANSION. 

In this Section, we show that our field theory is as thorough as the standard methods in statistical physics. In 
order to do this we compare our expansion with the standard Mayer expansion. To simplify the discussion, we first 
consider the case where the external potential is zero. 

To elaborate this comparison, we take the standard expansion of the grand potential in Mayer functions : /(r) = 
exp[— /3u(r)] — 1 and activity p z (28j. We further expand the exponential in terms of — (3v and obtain one and the sum 
of graphs with n > 1 lines in parallel representing the potentials multiplied by a factor 1/n!. This corresponds to an 
expansion in terms of the single potential introduced in Section [Vl with all possible connected unlabelled graphs with 
vertices corresponding to the activity. In the following, this expansion will be referred to as Mayer expansion. The 
Mayer and Feynman graphs have the same topological elements. They both include all possible connected graphs 
made of lines and points. To then state the equivalence between the two expansions, we must discuss the following. 
Firstly, although they have a similar topology, lines and vertices are associated with different quantities. In FT, 
points are n T -vertices and one body external weights, whereas in the statistical mechanics they represent the activity. 
Secondly, we need to compare the combinatory coefficients for the two expansions. Hereafter, we do not discuss the 
powers of a as in the final the graph is proportional to a and focus only on p z . 

First we discuss the powers in p z a 3 . From the previous section, Feynman diagrams are based on the n T -vertex 
functions associated with l/(p z a 3 ) n ^ 1 . On these vertices, interaction potential lines are attached which from relation 
(|39|) each contribute with (p z a 3 ) 2 (one instance is explicit and one comes from the definition of v). This factor can be 
distributed on the two vertices to which any line is attached. By doing this we associate a single power p z a 3 to each 
vertex and none to the lines. The role of the one body vertex has to be treated separately. In one case, the external 
weight v is attached to the ideal system propagator. From eq. ([37)) , the only factor p z a 3 of the ideal propagator is 
already distributed to the vertex inside the diagram. One can verify that there remains one factor p z a 3 associated 
with the external weight v and this term corresponds to the activity which should be at the end of an external line 
in the Mayer expansion. In the second case, the external weight is attached to a single interaction potential. This 
corresponds to two potentials in series and we can use the two body vertex. We can verify that here also we have the 
correct number of factors p z a 3 once they are redistributed on each vertex and that we retrieve the standard Mayer 
graph result. The final statement is that although in the Feynman expansion factors p z a 3 are associated both with 
lines and vertices, they can formally be redistributed in order to associate a single instance of this coefficient to each 
vertex. This corresponds to the Mayer diagrams expansion where the activity is associated to the points. 

The second aspect is that the n T - vertex functions are associated with the standard 1/n! factor. This is exactly 
the correct combinatory so as to obtain the non labelled graphs of the Mayer expansion, with identical rules for 
the symmetry of the graphs. One only needs to treat separately the case of the external weight, which includes the 
interaction potential when it is attached to the ideal propagator. It corresponds in the Mayer expansion to a single 
potential pending from a graph, the topological equivalent of an external branch in the Feynman expansion. The 
expected combinatory is found in this case too. 

The sum of all these results shows that the Mayer and the Feynman expansion are finally identical. The result can 
be extended to the system in the presence of an external field, indeed we have shown that any vertex function can be 
decorated by a factor e~ v . We have seen above that each vertex function can be associated to a factor p z a 3 , the 
multiplication of this factor by the exponential corresponds, in the liquid state physics, to the generalization of the 
activity in the presence of an external field denoted z* = p z e~ v, " mt in [28| . 

The foremost result of this paper shows that given the renormalization introduced in Section HV El the result for the 
diagrammatic expansion is simple and leads to the equivalence of the Feynman and Mayer graph expansions. We thus 
fulfil our main objective which is to define a FT capable of describing the system at a microscopic level introducing a 
simple and intuitive Hamiltonian. This confirms previous results where we have shown that our formalism reproduces 
two exact results which are the virial theorem [2!| and the contact theorem [30[ which will be discussed in more detail 
in the next Section. Note that in our formalism, there is no reference to Gibbs ensembles. The difference in number of 
degrees of freedom associated with the field description and the lattice spacing a calls for the renormalisation which 
we have introduced in order to reproduce the correct combinatorics for particles. 
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VII. DISCUSSION. 



From standard text books [27], |3l| , we know that the so-called classical statistical mechanics contains two basic 
properties governed by quantum physics. Namely, the thermal de Broglie wavelength, A, and the indiscernibility of 
particles which originates from N distinct particles a TV! coefficient in the partition function. These elements are not 
related to the interaction potential. In the present paper, we have shown that a simple local functional together with 
a renormalization procedure can reproduce these two properties. This procedure is not modified when an interaction 
pair potential is introduced in the Hamiltonian and consequently we can then demonstrate that the theory is equivalent 
to the usual statistical mechanics. We have shown that the local functional leads, in perturbation theory, to a simple 
combinatory of the fields. In each monomial term, the n fields are equivalent and their permutation is associated with 
the coefficient 1/nl. In other words, the local functional transposes to the FT the indiscernibility of particles. 

One characteristic of our FT is that we have been able to introduce a renormalization procedure through which 
all the results are finite and independent from, arbitrary lattice spacing, although there exists an infinity of coupling 
constants. Due to renormalization, the expression exp{— (3H [p]} is formal and we must consider that this quantity is 
defined by its series expansion around the activity and that some terms in this expansion are cancelled by counter 
terms; these are independent of the interaction potential, showing that they have no physical meaning but are 
originated only by a mathematical procedure. 

Achieving a microscopically faithful description shows that a simple FT is not necessarily associated with a coarse 
graining and can have a level of description equivalent to that of the standard statistical mechanics, in contrast to the 
common conceptions of this type of approach (32J ■ Indeed, the measure we have used does not require the introduction 
of any normalization constant in the partition function, necessary in the case of a coarse grained approach. 
We can also compare this FT with other microscopically exact field theoretical descriptions. Considering a field 
approach without using as a starting point the standard partition function, we deal with a renormalization that does 
not exist for field theories based on the Hubbard Stratonovich transform. On the other hand, our field is extremely 
simple and has an obvious physical meaning. This contrasts with the Hubbard Stratonovich type approaches, where 
we have to work in a complex plane with an auxiliary field for which it is rather difficult to introduce appropriate 
physical approximations. 

We also emphasize that FT is distinct from the DFT. Both approaches are based on the existence of a functional 
of the density. However, in the two formalisms, the correlations are treated in different ways [321. |33|. In the DFT, 
the form of the functional includes all correlations and fluctuations and we know that this functional exists but we 
ignore its exact form. Minimizing this functional yields the equilibrium density distribution. In contrast, in FT the 
functional is known and simple. The core of the FT formalism is to gradually account for the fluctuations when 
calculating quantities for the system in a perturbative expansion. One part of H[p] is F[p] which is formally like the 
free energy of the ideal system. A similar term Foft[< P >] is introduced in DFT, however it is important to point 
out the differences between F[p] and Fdft[< P >]■ F[p] is a functional of a field i.e. a fluctuating quantity where as, 
at the minimum, Fdft[< p >] is a function of the mean value of the fluid density. Moreover, we have mentionned 
earlier that exp{— f3H [p]} is essentially a formal expression. This illustrates one specificity of the FT : the fluctuations 
of the ideal term which basically represent the entropy must be considered on the same footing as the fluctuations 
related to the interaction pair potential. 

In this respect, although our FT is equivalent to standard statistical mechanics, the two approaches focus on different 
aspects of the correlations. This is the case when comparing standard approaches, but the discussion will be extended 
below for the case of our FT. We are convinced that having at disposal distinct formulations for a given quantity is 
indeed useful, possibly for acquiring a broader understanding. 



Hereafter we illustrate on three examples how FT leads to a new point of view on traditional quantities. 
In liquid state theory there are three classical expressions of the chemical potential. One of them corresponds to 



A. Examples 



34] 



ln(p(r)A 3 ) + In < exp(/3u(r)) > +V ext {r) = (ip. 



(41) 



A second traditional expression is given by [H, HH, [3^ | 



ln(p(r)A 3 ) - c (1) (T, [p];r) + V ext (r) = (3p 



(42) 
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where c^(T, [p]; i) is the single-particle direct correlation function 28, 3^,H(|. Finally, we also have a relation based 
on a charging process of the interaction potential [27| 

ln(p(r)A 3 )+p(r)^ d? J dr>(|r - r'|) 5 ( 2 )(|r - r'|;0 + V ext (r) = (3p (43) 

where ff^ is the pair distribution function [28] as a function of the charging parameter £. We note that all 

these expressions emphasize properties related to the potential, whether calculating the correlations of a quantity 
involving the interaction, or calculating the single-particle direct correlation function or alternatively considering a 
charging process of the interaction. 

The field theoretical description leads to a new expression which can be obtained by writing that the field is a 
dummy variable in the functional integral. This leads to the so called Dyson relations [3ll . and we obtain 

V/a 3 

< ln(p(i)A 3 ) > + 0v(i,j) < p(j) >a 3 + V &xt (i) = Pp. (44) 

Here the term related to the interactions is rather simple, it expresses the mean potential at a given point without 
taking into account the correlations. All correlations and fluctuations appear in the calculation of the average of the 
logarithm of the density field. This contrasts with a simple term like the logarithm of the average density, which 
appears in standard statistical mechanical expressions or in the DFT. As a consequence, differences in the description 
and a different organisation of the perturbation expansion in the FT, suggest that one should be able to elaborate 
new approximations. 

Let us consider now the so-called contact theorem which establishes an exact relation between the pressure p existing 
in a bulk phase and the value of the density profile p(0) at the wall enclosing the bulk material. This corresponds to 

Pp = p(0). (45) 

In so far as this relation is concerned, discussing the derivation of this theorem is an opportunity to emphasize 
the conceptual differences between the various approaches. We mention the kinetic theory of gases, in which this 
relation is the consequence of the mechanical equilibrium at the interface. In the case of DFT, the derivation is 
straightforward as we only need to write a displacement of the external potential : the interface, in two different 
ways. Another derivation [381 ] is obtained by integrating the BGY equations. In this case, a subtle integration of the 
correlations through the interface leads to the relation. Within our field theoretical framework, the key element is 
the local functional which is essential at different levels. It is crucial to obtain the density contact value present in 
the contact theorem [30| . but it is also necessary to cancel supplementary terms which appear in the demonstration. 
In this respect, specific relations of the field theory are also required, namely the Dyson type relations [371 ] . 



Now, in a third example, we illustrate one of the main aspects of the FT, i.e. the existence of an intricate coupling 
between counting (entropy) and interaction. Let us study the interfacial properties of ionic fluids. From the point of 
view of the interactions, we know that the important quantity is the charge, the difference of densities of each species. 
However, this system can also be viewed as a peculiar mixture which has a specific condition due to electroneutrality. 
From this point of view, we have two terms in the ideal functional describing the indiscernibility for each ion. Thus 
the natural fields are the densities describing each ion. In [3^, H(|, we show, in the specific instance where the natural 
fields for the ideal and for the interaction term are distinct, that the perturbation theory leads to a coupling of the 
charge and of the total density field due to the local ideal functional. This has direct consequences. 
For the simple neutral interface, we show that there exists a depletion for the quadratic fluctuations of the charge. 
Then, the entropic coupling between the charge field and the total density field predicts a non trivial profile on the 
total density [H, • We can verify that the contact value of this total density profile satisfies the exact condition of 
the contact theorem, for the pressure calculated at the same level of approximation. We have used this phenomenon 
to analyse the anomalous behaviour of the differential cap acitance as a function of the temperature [ill, E3] which 
has been thoroughly discussed recently in experiments |43l . l44j , numerical simulation (4f| and theoretical approaches 
HE S3, 0) El] • The interest of our analysis is that it provides a simple interpretation and understanding for this 
phenomenon, associating the decrease in the capacitance with the depletion of the ionic density at the interface at 
low temperature and providing the physical origin of this depletion. 

Moreover, the more detailed account of these entropic effects is fundamental in the case of asymmetric, in valence, 



electrolytes. In [42j, we have tested our FT by comparing with the results of numerical simulations [50| and shown 



that the theory accounts for all main qualitative properties of the phenomenon, in comparison with other approaches 
51] which although currently more quantitative fail to take certain features into consideration. 
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VIII. CONCLUSION 



In this paper, we present a field theory describing classical fluids at equilibrium at the same level as the standard 
statistical mechanics. We introduce a real physical field and construct the Hamiltonian in the spirit of the QFT. 
This functional includes interactions and a local functional representing the ideal system. The latter characterises our 
approach and has been thoroughly discussed. In particular, we show that it provides, for the FT, essential ingredients 
in relation to quantum mechanics. The equivalence of our theory with standard statistical mechanics is shown by 
establishing that the Feynman expansion of the FT is equivalent to the standard Mayer expansion. The approach is 
original in that it is not a simple mapping of the standard partition function like other field theories. Consequently, 
it requires a renormalization which we describe. Its basic interest is that the theory remains simple and intuitive like 
phcnomenological field theories. 

Establishing a field theory which is both a simple and exact representation of the statistical mechanics has many 
advantages. We present possible applications. 

Some are related to the FT formalism. We can, for instance, use powerful tools such as discussions in terms of 
symmetries of the system, of the fields Also, the fact of having a field variable at the microscopic level should 
allow for natural bridging with the mesoscopic intuitive approaches which also adopt the field theory description. An 
example can be found in [53| where a mesoscopic Hamiltonian is presented. 

Another aspect is that this formalism treats fluctuations in a different way. This type of approach would help 
elabora ting small systems, where fluctuations can have the same magnitude as the quantities characterising the 
system |32j |. 

Finally, we have also shown that there is an emphasis on correlations associated with entropic effects. Such 
emphasis should shed new light on the description of ionic systems, or mixtures. For instance, we believe that 
the emphasis on the correlations between charge and total density could add to the understanding of criticality in 
ionic systems. For such systems, the potential couples the charge, whereas criticality characterizes a phenomenon 
on the total density. Another system of interest in the field of the double layer is the study of asymmetric 
in charge electrolytes, which exhibit polarization phenomena even in the vicinity of neutral interfaces. As op- 
posed to asymmetric in size ions, this phenomenon is not intuitive. The difference of density of anions and cations 
for these asymmetric systems seems to be the origin of such phenomena as a consequence again of entropic effects 54 1 . 
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APPENDIX A: BEYOND THE IDEAL SYSTEM SADDLE POINT 



Beyond the saddle point, we can compute the integral eq. ||5J) taking into account the fluctuations of the field, on 
each lattice site we expand the density field as p = p + 5p, in this case the logarithm of the partition function is 



V 

In 5q = pV H — t In 



die~* 2/2 exp ( ^ 



(-l) n+1 t" 



\ n — 3 



n(n - l)(/5a 3 ) n / 2 



(Al) 



where t = (p — p)a 3 / ' \J pa 3 . Expanding the last exponent, which makes sense in the limit of large pa 3 , we find we 
have to calculate Gaussian integrals: t 2n e~ f l 2 dt — \f2ii{2n — 1)!!. The result can be written 



V O 

In H =pV+ ^{pa 3 } = pV 



pa 3 



ip[pa 



(A2) 



with 



ip[pa- 



1 

-ln(27r/5a 3 ) + ^ 



L ^ 2 (pa 3 )^" 1 ) 



i i 

s~i/~f 3 \ 2 



161 1 



(pa 3 ) 2 5760 (pa 3 ) 3 



(A3) 
(A4) 



with the exclusion of the first term, tp is a power series of 1/pa which is asymptotically convergent for large pa , for 
which the first values of the coefficients dz, are given on the second line. 
The expression in the presence of an external potential V ext (ri) is 



tp[pe 



-V*(rO ft 3l 



1 1 °° 

= ln(2^a 3 )--V- t (r l ) + £ 



L=1 



(pa 3 )^" 1 ) 



(A5) 
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